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ABSTRACT 

We show how the cosmological constant can be estimated from redshift surveys at 
different redshifts, using maximum- likelihood techniques. The apparent redshift-space 
clustering on large scales ( ^ 20 /i _1 Mpc ) are affected in the radial direction by infall, 
and curvature influences the apparent correlations in the transverse direction. The 
relative strengths of the two effects will strongly vary with redshift. Using a simple 
idealized survey geometry, we compute the smoothed correlation matrix of the redshift- 
space correlation function, and the Fisher matrix for Qm and Oa- These represent 
the best possible measurement of these parameters given the geometry. We find that 
the likelihood contours are turning, according to the behavior of the angular-diameter 
distance relation. The clustering measures from redshift surveys at intermediate-to-high 
redshifts can provide a surprisingly tight constraint on We also estimate confidence 
contours for real survey geometries, using the SDSS LRG and QSO surveys as specific 
examples. We believe that this method will become a practical tool to constrain the 
nature of the dark energy. 

Subject headings: cosmology: theory — galaxy clustering — large-scale structure of 
universe 



1. Introduction 

The dark energy, which is the cosmological constant A in its simplest form, currently have 
turned out to possibly be a dominant component in the universe (e.g., Bahcall et al. 1999). It is 
one of the central issues in cosmology to reveal the quantitative nature of this mysterious form of 
the energy. One of the mysteries of the cosmological constant is its smallness. There are no evidence 
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detected on the Earth for existence of the cosmological constant. The cosmological constant is so 
small that it only affects the phenomena on cosmologically large scales. 

There are several traditional tests of the cosmological constant. The expected frequency of 
multiple image lensing events for high-redshift sources is quite sensitive to the cosmological con- 
stant (Fukugita, Futamase, & Kasai 1990; Turner 1990). Luminosity-volume, and redshift-volume 
relations can also be used to measure the geometry of the universe to constrain the cosmological 
constant (Rowan-Robinson 1968; Loh 1988). Alcock & Paczyhski (1979) proposed an evolution- 
free test for the cosmological constant using statistically spherical objects. The type la supernova 
Hubble diagram is used to constrain the mass density parameter and the cosmological constant 
(Schmidt et al. 1998; Perlmutter et al. 1999) through the luminosity distance di- Acoustic peaks of 
Cosmic Microwave Background (CMB) anisotropies constrain the curvature of the universe (e.g., 
Hu, Sugiyama & Silk 1997). Recent observational developments of the type la supernova and the 
CMB anisotropies (Balbi et al. 2000; de Bernardis et al. 2000) suggest a flat, low-density universe 
with positive cosmological constant, f^M ~ 0.3, SIa ~ 0.7. 

The effect of a cosmological constant on the clustering properties of objects in the nearby 
universe (z <C 1) is so weak, that redshift surveys have not been used to constrain Oa so far. 
As the depth and the sampling rate of redshift surveys increase, redshift-space clustering depends 
on the cosmological constant through the cosmological redshift distortions (Ballinger, Peacock & 
Heavens 1996; Matsubara & Suto 1996; Matsubara 2000). Several applications of this effect are 
proposed (Nair 1999; Nakamura, Matsubara & Suto 1998; Popowski et al. 1998; Yamamoto &; 
Suto 1999; Yamamoto, Nishioka & Suto 1999). To maximally extract the cosmological information 
from the survey data, the likelihood analysis combined with a data reduction technique like the 
Karhunen-Loeve transform has been quite successful at low redshifts (Vogeley & Szalay 1996; Szalay, 
Matsubara & Landy 1998; Matsubara, Szalay & Landy 2000). We expect it to be just as useful at 
intermediate-to-high redshifts here. 

In this Letter, we combine the two methods, i.e., the likelihood analysis on pixelized data 
and the cosmological distortions. Specifically, we compute the Fisher matrix for simple geometries 
to illustrate how the cosmological distortions constrain the cosmological constant and the density 
parameter in a given redshift survey data. 

2. From Correlations to Fisher Matrix 

To generically investigate how a given redshift survey can constrain the cosmological constant, 
we construct a rectangular box, in which objects like galaxies or quasars in redshift space are 
observed. We will use smooth pixels rather than hard cells in order to make our calculations 
numerically more efficient. We apply a Gaussian smoothing window to the objects in the survey, to 
get a smoothed estimate of the local density. With a sufficiently large smoothing radius, we do not 
have to deal with the nonlinearities of the density field. The Gaussian smoothed cells are placed 
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on lattice sites in the box. In this way, the smoothed density vector pi on discrete lattice sites 
labeled by i is considered as our fundamental data to be analyzed. In the following we assume the 
mean value of the density vector (pi) is known so that we can define a density-fluctuation vector 
<k = Pi/{pi) - 1. 

In standard theories of structure formation, the linear density field is a random Gaussian pro- 
cess. In this case, all clustering properties of the universe are represented by two-point correlations. 
Thus, a correlation matrix 

Qj = {didj) , (1) 

theoretically specifies all the statistical information for a given data set. This matrix is related to 
a smoothed two-point correlation function plus a shot noise term: 

Ctj = j d 3 Sl d 3 s 2 W(s t - Sl )W{s 3 - s 2 ) s 2 ) + J d 3 s W{s t - s)W{s 3 - s)/n(s), (2) 



where W(s) = exp[— s 2 /(2R 2 )]/(\^2ttR) is a Gaussian smoothing window, £ is a two-point correla- 
tion function, and Sj is the position vector of a lattice site i, and n(s) is the mean number density 
field. The position vectors are all in observable redshift space which consists of the redshift z and 
the angular position (6,<j>). In the distant-observer approximation, we can approximately use the 
Cartesian coordinates in redshift space and the two-point correlation function is a function of the 
relative vector s\ — s 2 . The correlation function does depend on the direction of this relative vector, 
because of the redshift distortions. An analytic form of the linear two-point correlation function 
in redshift space including high-redshift effects is given by Matsubara & Suto (1996). One of the 
equivalent forms given in Matsubara & Suto (1996) is (see also Ballinger, Peacock & Heavens 1996) 

- s 2 ) = b\z)D\z) j -0^e ik < x ^ [l + [3(z)k 3 2 /k 2 ] 2 P(k), (3) 

where P(k) is the linear mass power spectrum at z = 0, D{z) is the linear growth rate normalized as 
D(0) = 1, b(z) is the bias factor at redshift z, and (5{z) is the redshift distortion parameter, which is 
approximately related to the redshift-dependent mass density parameter as (3{z) ~ £lM 0M (z)/b(z). 
In the above equation, the third axis is taken as the direction of the line of sight. The vectors x\ 
and x 2 are the comoving positions of the two points which are labeled by s\ and s 2 in redshift space 
(see Matsubara &; Suto 1996), i.e., they are related by a comoving distance-redshift relation and 
the spatial curvature of the universe. Assuming the distant-observer approximation, a Gaussian 
window function, and that the mean number density is effectively constant, n, the equation (2) 
finally reduces to 

dj = -b 2 (z)D 2 (z) [°° ^P(k) r sinfldflexp {-k 2 R 2 [(c,, 2 - c ± 2 ) cos 2 9 + c ± 2 ] } 
2 Jo 2tt J 

x (l + ficos 2 6) Jo(kxsm6 x sin9)cos(kxcos9 x cos9) 

exp [-x 2 /(AR 2 )] 

^/ 2 {2Rfn 1 ' 
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In this equation, c\\(z) = Hq/H(z), c±(z) = Hqsk(z)/ z are the distortion factor parallel and 
perpendicular to the line of sight, respectively, where H{z) and sk(z) are the Hubble parameter 
and the comoving angular diameter distance at z, and Hq is the Hubble's constant. A line-of-sight 
component of the redshift-space distance sm between the centers of i-cell and j-cell is related to 
that of the comoving distance xy by xy = C||(z)s||. Similarly, for a component perpendicular to 
the line-of-sight we have x± = c±(z)s±. In this notation, the quantities in equation (4) can be 
written as x = (c|| 2 S|| 2 + c± 2 s± 2 ) 1 ^ 2 , and 9 X = cos _1 (c||S|i/x). Integration over 6 remains because 
a spherical Gaussian smoothing kernel in observable redshift space is no longer spherical but is 
ellipsoidal in comoving space. The second term in equation (4) is a shot noise term, convolved with 
the Gaussian kernel. 

Once the correlation matrix can be theoretically calculated in any cosmological model, one 
can calculate the Cramer-Rao bound which gives an estimate how well the model parameters can 
be measured. This is one of the most powerful results in estimation theory (Therrien 1992). The 
Fisher information matrix is a key quantity in this theory (Kendall & Stuart 1969): 

where L(d; 6) is a probability distribution for the data vector d, which depends on a vector of model 
parameters 6. In our case, the data vector is density fluctuations on lattice sites and the model 
parameters are the cosmological parameters. The Cramer- Rao bound states that the maximal 
likelihood estimate constrains the model parameters with a minimum variance 

(e a e?)>{F- l ) a ^ (6) 

where F~ l is the inverse matrix of F. When the number of data, i.e., the dimension of the data 
vector is very large, the Cramer-Rao bound (6) becomes equality. A contour [0 a F a p6p\ 1 / 2 = A in 
parameter space gives a concentration ellipsoid, which indicate regions where the likelihood density 
for model parameters are most concentrated. The threshold A = 1,2,3 corresponds to maximally 
attainable confidence levels of la, 2a, 3a, respectively, in a likelihood analysis, if the likelihood 
function is Gaussian. The concentration ellipsoids are useful even when the likelihood function is 
not Gaussian to give a rough idea of the spread of the density (Therrien 1992). 

As we are interested in a linear density field, the probability distribution of the density field is 
considered to be Gaussian so that the likelihood function has the form, 

-21nL = In detC + d T C~ 1 d + const. (7) 

where C is the correlation matrix which depends on model parameters, and d is the data vector. 
In this case, the Fisher information matrix reduces to (see, i.e., Vogeley &; Szalay 1996) 

F aP = ^Tr(C- 1 C, a C- 1 C,p), (8) 

where C jQ , = dC/d9 a , etc. Thus, the Fisher matrix or the concentration ellipsoids for any model 
parameters are straightforward to calculate from the correlation matrix of equation (4) . 
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3. Results for a Simple Cubic Box 

In this Letter, the mass density parameter Qm and the normalized cosmological constant U\, 
both at the present time, are the model parameters to be constrained. These parameters has the 
primary importance in high-redshift clustering distortions. For simplicity, the power spectrum P(k) 
and the bias parameter b(z), on which the correlation matrix of equation (4) also depends, are fixed 
throughout. We use the cold dark matter-type spectrum with a fixed shape parameter T = 0.2 and 
a fixed normalization as = 1. The growth factor D{z) and the distortion parameter (3{z) are the 
functions of £Im and Throughout this Letter, we take a fiducial models (Qmj^a) = (0.3,0.7) 
at which the Fisher matrix is evaluated. 

Although the most natural choice of the length unit system would be comoving coordinate 
system, we should not use this system, because they are not actually observable in redshift surveys, 
and depend on the cosmological models which we are seeking. Thus, we use the coordinate system 
like (z, 9, 4>) in polar coordinates. The clustering scale we are interested in is too small in figures 
when the distance is represented by z itself, so that we invent a new radial coordinate s = cz /Hq ~ 
2997.9 z, i.e., the linear extrapolation of distance-redshift relation for z <C 1. For example, a 
redshift interval Az = 0.1 around any redshift z corresponds to As = 300. Although the unit of 
this coordinate system is still h~ 1 Mpc , we use a new notation /i _1 Mpc 2 to avoid a confusion with 
the comoving coordinate system. 

In this coordinate system, a 200/i _1 Mpc z cubic box is considered to obtain generic estimates 
for Cramer-Rao bound, and we compute the Fisher matrix for this sample, varying the mean 
redshift z of this box. The density fluctuations are sampled on regular 10 x 10 x 10 lattice sites in 
the box. The Gaussian smoothing radius is set as R = 10/i _1 Mpc z . 

Figure 1 shows the concentration ellipses of the 200/i~ 1 Mpc 2 box placed at redshift z = to 
6.0. The contour lines correspond to A = 1, 2, 3 which indicate the maximally attainable confidence 
levels of la, 2a, 3a when one performs a likelihood analysis for these samples as we described in 
the previous section. The shot noise is neglected in this Figure, while it would be difficult to 
reduce shot noise for z > 2 in reality. As is well known, a low-redshift sample (z ~ 0) only 
constrains the mass density parameter through its dependence of the redshift distortion parameter 
(3 ~ Q^/b. Increasing the mean redshift, the concentration ellipses rotates clockwise and the 
major axis becomes shorter, and thus the cosmological constant becomes constrained. The higher 
the mean redshift is, the more the cosmological constant is constrained. Around z ~ 1.7, the 
concentration ellipses begin to rotate counterclockwise. This is consistent with the fact that the 
angular diameter distance-redshift relation turns over at z ~ 1.7 in our fiducial model (^m, ^a) = 
(0.3,0.7). 

On one hand, the number density one can sample is smaller for high-redshift objects, which 
dilute the constraints on cosmological parameters. On the other hand, there is a larger volume to 
be sampled than for low-redshift objects. To obtain the concentration ellipses in realistic samples, 
one should take into account both the shot noise effect and the total volume in a given sample. We 
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Fig. 1. — Concentration ellipses from the Fisher matrix for generic boxes. Contour lines correspond 
to confidence levels la, 2a, 3a attainable from a single 200ft,~ 1 Mpc z box placed at redshift z = 
to 6.0 as indicated in each panels. This figure contains no shot noise, and a bias factor of 1. 
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again set the 200/i~ 1 Mpc 2 boxes and estimate the Cramer-Rao bound with the shot noise effect 
included. In Figure 2, the bound for the normalized cosmological constant S7a is plotted. In the 
left panel, the volume is given by just one 200/i -1 Mpc 2 box as in the case of Figure 1. The shot 
noise is varied as (20 /i~ 1 Mpc z ) 3 n = 0.1,0.3, 1,3, 10, and oo, from top to bottom lines. In the right 
panel, the Cramer-Rao bound is scaled by the number of independent 200/i -1 Mpc z boxes in a 7r 
steradian region with a redshift interval z/2 around z, to obtain a rough idea of how our error 
bound is affected by the survey volume. We can see how densely the objects should be sampled to 
constrain the cosmological constant with a certain accuracy both for a sample with a fixed volume 
and a sample with a fixed solid angle. 



4. Possibilities for Realistic Survey Volumes 

We have considered several different survey layouts for both galaxies and quasars. The best 
survey to perform these tests seems to be the Luminous Red Galaxy (LRG) sample of the Sloan 
Digital Sky Survey (SDSS). This sample consists of 100,000 galaxies selected for spectroscopic 
observations on the basis of their very red rest frame colors, using photometric rcdshifts, down to 
a limiting magnitude of r' = 19.5. They form an approximately volume limited sample, where the 
outer edge lies at around z = 0.45, and the total surface area is 10,000 square degrees. 

We consider this geometry as a composite of the generic 200/i _1 Mpc z boxes at the mean 
redshift z = 0.3. There are about 220 boxes out to z = 0.45 in a ir steradian region, so that the 
shot noise is approximately given by (20 h~ 1 Mpc z ) 3 n = 0.5. We assume two bias factors, 6 = 1.5 
and 6 = 2. The resulting concentration ellipses are shown in Figure 3. We can see the result is quite 
remarkable. The Cramer- Rao bound for the cosmological constant is only [(-F _1 )aa] 1 ^ 2 = 0.04 for 
6 = 1.5, and [(F _1 )aa] 1 / 2 = 0.03 for 6 = 2. This shows that the shot noise level and the depth 
of the survey volume are suitably balanced to constrain the geometry of the universe in the SDSS 
LRG survey. 

We have considered various QSO surveys to possibly measure the dark energy at higher red- 
shifts. Unfortunately, the currently ongoing QSO redshift surveys, like Sloan Digital Sky Survey 
(SDSS, i.e., York et al. 2000) and 2dF QSO redshift survey (2QZ, i.e., Boyle et al. 2001), have lower 
sampling rates for QSOs, n ~ 10" 3 /(40/i" 1 MpcJ 3 . They typically give the Cramer-Rao bound of 
order A^a ~ 1, almost regardless of the smoothing radius. To constrain the cosmological constant 
with QSO surveys, one should sample QSOs more densely than these current QSO surveys. This 
fact is in agreement with Popowski et al. (1998) who analyzed nonlinear clustering to constrain the 
geometry of the universe and indicated an advantage of a dense sampling. 



-8- 




Fig. 2. — Cramer-Rao bound for the normalized cosmological constant J7a as a function of the 
mean redshift z. Shot noise is varied as (20 /t _1 Mpc z ) 3 n = 0.1,0.3,1,3,10, and oo, from top to 
bottom lines. Left panel: a single 200/i~ 1 Mpc 2 box. Right panel: For a survey of 7r steradians, a 
redshift interval of z/2, centered at z. The number of independent boxes is increased by the volume 
of the shell, correspondingly the accuracy improves. A bias factor of 1 has been used throughout 
this figure. 




Fig. 3. — Concentration ellipses corresponding to la, 2a, 3a confidence levels for approximate 
geometries of the 100,000 galaxies in the SDSS LRG (Luminous Red Galaxy) sample. Dotted lines 
assume a bias factor of b = 1.5, solid lines has 6 = 2. 
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5. Discussion 

We have shown that large-scale clustering of galaxies at intermediate redshifts z ~ 0.5 is 
surprisingly suitable for constraining the cosmological parameters of Qm and f^A, and thus the 
geometry of the universe. The QSOs in currently ongoing surveys are too sparse to give comparable 
constraints. 

The apparent redshift-space clustering method used in this Letter is a completely self-contained 
test for f^M and Oa- The results from this method can be further combined with any of the other 
independent tests to obtain stricter constraints, or to check a consistency of our standard picture 
of the cosmology. 

One can use Figure 2 to aid designs of future surveys at various redshifts. The lines indicate 
the statistical uncertainty in the cosmological constant corresponding to different sampling rates. 
They should be scaled, noting that the Cramer-Rao error bound roughly scales as the inverse of 
the square of survey volumes. 

In this work, we have only considered two parameters S7m and S7a- We still need to measure 
the evolution of bias parameter, which is not obvious. Moreover, there is a possiblity that the dark 
energy has a more complex behaviour than the cosmological constant (Wang et al. 2000). There 
are many other cosmological parameters, like baryonic density f2t>> primordial spectral index n, the 
neutrino mass density VL U , etc., which more or less depend on the apparent redshift-space clustering. 

AS acknowledges support from grants NSF AST-9802 980 and NASA LTSA NAG-53503. We 
would like to acknowledge useful discussions with Dan VandenBerk, Daniel Eisenstein and Adrian 
Pope. 
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